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^ Abstract 

The luminosity measurement at the projected International Linear e'^e Collider ILC is planned 
^ ■ to be performed with forward Bhabha scattering with an accuracy of the order of 10~'^. A theo- 
■ retical prediction of the differential cross-section has to include one-loop weak corrections, with 
leading higher order terms, and the complete two-loop QED corrections. Here, we present the 
weak part and the virtual one-loop photonic corrections. For the photonic corrections, the ex- 
pansions in e = (4 — (i)/2 are derived with inclusion of the terms of order e in order to match 
the two-loop accuracy. For the photonic box master integral in d dimensions we compare several 
different methods of evaluation. 
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1 Introduction 



Bhabha scattering 

e"(pi) + e+(p4) e~{-p2) + e^{-p^) (1-1) 

was one of the first processes calculated in quantum theory [1]. The complete virtual electroweak 
one-loop corrections have been first calculated in [2], later also in [3-11]. By now, Bhabha 
scattering may also be calculated with automated tools for the evaluation of Feynman diagrams 
and cross-sections as e.g. Feynarts [12, 13], grace [14] and alTALC [15]. The electroweak 
corrections have to be considered together with hard bremsstrahlung corrections, which usually 
are calculated by Monte Carlo programs; see [16-20] and references therein. Dedicated studies 
for experimentation at LEP may be found in [21, 22] and references therein. 

The preparation of the e+e^ linear collider project ILC (formerly also TESLA [23], and cor- 
responding projects of other regions) triggered again some interest in both wide angle and small 
angle Bhabha scattering. The latter might allow to determine the luminosity with an unprece- 
dented accuracy of 10^^. For this, one needs theoretical predictions beyond one-loop accuracy 
in the extreme forward scattering region where the cross-section peaks due to the kinematical 
singularity of the photon propagators, while the pure weak corrections might be sufficient in 
one-loop approximation (with leading higher order terms a la [7]). If a so-called Giga-Z option 
will be realized, high Bhabha event rates are to be expected in the Z resonance region also for 
larger scattering angles. 

We write the matrix element squared 

+ + 23fJ(A<(°)*A<(2)) ^ 23fJ(A<f + (1.2) 
V ' 

where A^(*) is the contribution to the i-loop order and the subscripts 7 and 77 indicate the emis- 
sion of one or two photons. 

The QED contributions dominate by far and two-loop corrections are also needed. Sev- 
eral projects to determine them in a systematic way are underway (see [24-29] and references 
therein). A program with two-loop accuracy has to include also the complete one-loop matrix 
elements squared, often regulated by an expansion in e = (d — 4)/2, with a careful treatment 
of the resulting finite terms in e.' For this, one may express the Feynman diagrams by scalar 
master integrals, which then have to be known up to some positive order in e. Thus, one has to 
go beyond the usual technical demands of a pure one-loop calculation. 

In this article, we give a concise description of our approach to the one-loop contributions 
for a two-loop calculation of massive Bhabha scattering. Introductory, we present in Section |2l 

'a similar program was performed in [30, 31]. 
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electroweak predictions which were obtained for the ILC study [23, 32, 33]. The expressions for 
the pure QED corrections up to order e in terms of a few scalar master integrals are derived in 
Section|3l Here we retain the exact dependences on the electron mass. The scalar master integrals 
are discussed in Section IH For the box master integral we compare several, quite different 
expressions, which are derived with the aid of a difference equation, a system of differential 
equations, and the Mellin-Barnes technique, respectively. We close with a short Summary. 

2 Electroweak one-loop corrections 

One-loop corrections are the virtual part of the 0{a^) terms in (11.21) . The calculation of elec- 
troweak corrections to Bhabha scattering with the automated tool a I TALC has been described on 
several occasions [15,32-35]. alTALC [15] uses the packages DIANA V.2.35/QGRAF 2 [36,37] 
for the creation of the one-loop matrix elements, FORM 3.1 [38] for their expressions in terms 
of scalar s, and LoopTools 2.1/FF [39,40] for the numerical evaluation, including also soft 
bremsstrahlung. In one respect we had to go beyond LoopTools 2.1: In order to evaluate 
cross-sections in the neighbourhood of the Z resonance peak, one has to use Breit-Wigner prop- 
agators in the s-channel, replacing m| by ml = m?z~'^^z^ z- Accordingly, the box function 
Do(t, s, mo) = Dq^iti^, m^, m?, m?, t, s, A^, m^, ml, m^), if used with (7, Z) in the s-channel, 
has been modified as in [41]: 



The diagram is shown in Figure IZTI 

We use the LoopTools conventions s = (pi+P4)^,t = (pi+P2)^,pf = = rn^-,^ = rn, 
with the following definitions for the x variables 





(2.1) 



Xs = —K{s + ie,mi,m4), 
X2 = -K{ml,mi,mo), 
Xs = -K{ml,mi,mo), 



(2.2) 
(2.3) 
(2.4) 



and the definition of the fC-function (with one of the arguments being complex) 




(2.5) 



—1 z = {m — m'Y- 
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Table 2.1: Differential cross-sections in pbarn for Bhabha scattering at y/s = mz- Bom contri- 
bution and the 0{a) correction are shown; the maximum soft-photon energy is -yi/10. 



The photon mass is A in LoopTools. The use of this scalar box function is not only necessary in 
order to regulate the 7Z box contribution, but also for a proper compensation of the correspond- 
ing soft photon infra-red divergencies which are proportional to the Born cross-section with a 
Breit-Wigner Z propagator. The implementation is done in the file fort ran/ src/dOwddO . F 
of the a I TALC package. 

In Tables 12.11 and 12.21 we provide numerical sample outputs at typical energies for several 
scattering angles. The input quantities as well as the treatment of soft photons are exactly the 
same as in [42]. The cross-section peak in the forward direction, due to the photon exchange in 
the t-channel. In this kinematic region, the pure photonic corrections will be dominating and we 
have to treat them with higher accuracy than the rest of the electroweak corrections. For the one- 
loop corrections, this means a determination of p as part of the 0{a^) terms in (11.21) . Here 
one needs the QED one-loop functions including terms of order e because their interference with 
other terms of order 1/e contributes to the finite cross-section. This will be the main concern of 
the rest of this article. 




Figure 2.1: Four-point function with complex mass mo and photon mass regulator A. 
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Table 2.2: Differential cross-sections in pbam for Bhabha scattering at ^ = 500. Bom con- 
tribution, the 0{a) correction, and also a QED prediction are shown; the maximum soft-photon 
energy is ^/s/10. 

3 The massive QED cross-section in d dimensions 

The ten diagrams of Fig. (13.11) are the one-loop contributions in pure QED. 
We decompose the full one-loop matrix element as follows: 



Ml = 


[it, ® 7m] ^1, 




M2 = 






Ms = 






M4 = 






M, = 


([^4 ® 1] + [1 ® ^2])F5, 




Me = 


([7m^4 ® 7/.] + [7/. ® 7m^2 


])i^6, 


M7 = 


[1 ® 1] Fj, 




Ms = 


[iplu ® 7^/7/.] ^8, 




M, = 


([7m7;.^4 ® 7;.7,J + [7^7!. 


® lult,i>2])F9 



The notation is short-hand for the s and t channel matrix elements: 

M^^ = 0,®OfFk{s,t) 

= v{p4)0iu{pi) ■ u{-p2)0fv{-ps)Fk{s,t), (3.2) 

= 0,^Op[-F,{t,s)] 

= v{pi)Oev{-ps) - u{-p2)0pu{pi) [-Fk{t, s)] . (3.3) 

Crossing the diagrams from the s-channel to the t-channel results in the exchange of s and t 
and in an overall sign change due to Fermi statistics; see Equation (13.31) . In general the first six 
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Figure 3.1: One-loop diagrams for the process e+e — e+e 



amplitudes are independent while F7 to Fg can be expressed in terms of them, but in slightly 
different ways for the various cases under consideration. For this reason we list them all. 

With the form factors Fk{s,t) one may determine the contributions from 2'Si{A4^^^ Ai^^^) 
and from |A^(^)p to the differential cross-section (11.21) . The interference of Ai^^^ and A4^^^ 
yields e.g.: 



da vra^ 



dcosO 2s 



J2[B,{s,t)F^{s,t) + B^{t,s)F^{t,s)] 



(3.4) 



with 



B2 



Bs 
S4 
B5 
B, 



s <! 4(1 + v)"" + 2{l-v + -)- Avz^ + (1 - ^)z^) 

V 2v 

(1 + v)[{l + vf + v^] + ^[(1 + vf + v{3 + 5v) - 
s 1 8(3 + 3t; + 5v^) + - + 8(3-5v + -)z^ + 2(5 - -)z^ 

I V V V 

I -4(1 + vy - (2 - + -)z^ + (3 - 4t; + hz" + (1 - -f-)^^ 

1 3 1 

ms { -2(<6v + 4t;2 --) + {2 + Sv- -)z^ - (2 - -)^^ 

V V V 

ms { -4(3 + -) + 2(1 + 2v)z^ - (2 - -)z^ 

V V 



Br = 


( 

■{ 


Bs = 


■{ 




ms 



and 



t 

s 

Aw? 



V ^ Av^ 2v 

V V V 

-8(3v + V - -) - 4(3 -8v + -)z^ - 2(4 - -)z^ \ , (3.5) 

V V V 



(3.6) 
(3.7) 



With the same formula (13.41) . the Bom cross-section, arising from | A^'^^^ | , is obtained with: 

1 



Ff-"(s,t) = -, (3.8) 



s 



Ff^^^'is^t) = for J > 1. (3.9) 

The contributions from [A^^^^p to the cross-section are rather lengthy and not shown here ex- 
plicitely; they will be provided on the webpage [43]. There we give also the expressions for the 
corresponding interferences in d dimensions. 

Before determining the form factors Fk{s, t), we discuss now the various contributions. As 
mentioned we may restrict ourselves to the s-channel diagrams Dl, D3, D6, D8, D9: 

F,{s,t) = {Ft'' + Fr + F^-) , j = l,...A (3.10) 

The self-energy contributes to Fi only: 

In a theory with several fermion flavors (with different masses mj), one has to sum this term 
over all flavors. The vertices contribute to Fi and F^: 

pvert ^ Ff^'^ + F^^, (3.12) 

with F/^^ = F^*^. The two form factors for A^s in (13.11) are also equal but contribute to different 
structures there. The situation for the box diagrams is a little more involved: 

Ff^' = C6(Ff8 + Ff9), (3.14) 

and the Cb will be given in (13.331) . As mentioned only six of the nine form factors are independent. 
For the direct box diagram D8 we find the following relations: 

= 4m2Ff8 + 2mF5^8^ 
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Fg^s = mF^^'-^F,^'. (3.15) 

There are further relations between the form factors Fj^^ from the direct box D8 and the Fj^^ of 
the crossed box D9: 

Ff« = -Ff9 + (4-6d)F3^^ 
^ ^2^9 - (4 - 2ci) ^ 

pD8 ^ pD9 

pDS ^ _ pD9 

Fg^s = -2dm F^^ - F^^ + d F^\ 

F^^ = F(f9-4mFi^^ 

pm ^ _ (4 - 2rf) F8^^ 

pD8 _ _ pD9 

Fg^S = Fg^^. (3.16) 

Inverting relations (13.161) . diagram D9 is obtained from D8 by exchanging t and u. As mentioned, 
the t-channel box D7 may be obtained from D8 by simply exchanging t and s (and an overall 
sign). Subsequently, diagram DIO results from D7 by inverting again (13.161) and exchanging now 
s and u. As a consistency check, one can additionally obtain DIO from D9 by s, t crossing. The 
inversion of the first six relations of (13.161) yields Ff^ to F^^: 

pD9 ^ _pD8^^^_Q^^pD8^ 



F2^y = F^"" - (4 - 2d) F, 
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TpD9 _ pD8 
^3 — -^3 ' 
pD9 ^ _pD8 

F^^ = -F^^-2dF^\ 

Fe^^ = -F6^8-4F9^^ (3.17) 

In a next step, one gets in combination with (13.151) : 

F,^9 = -Qdm^F^^-2mF^'' + 2dmF^\ 
F^' = m'F^\ 

Ff' = mF^'-^-F^\ (3.18) 

We see that the relations for F-j in terms of amplitudes Fi to Fg are different for diagrams D8 and 
D9. 

What remains now is to determine one form factor for the self-energy, two form factors of 
the vertex, and six form factors for one of the four box diagrams. This will be done in two steps. 
First, we collect the form factor contributions from the Feynman diagrams Dl to DIO, and in a 
second step we have to add up additional contributions Fj*'*^ arising from counter term insertions 
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into the one-loop diagrams. The latter are formally of higher order, but it is reasonable to discuss 
them here. So, effectively, (13.101) has to be replaced by 



^ / N 2e' 



■pself I rpvert , rpbox , / rpself,r . rpvert,r . rpbox,r 



(4vr)d/2 

j = l,...,9. (3.19) 

Additionally, charge renormalization Se/e will give an overall factor, and there are also contri- 
butions Ff, F^'^ from wave function renormalization. Both will be discussed in Section |3]2| 

3.1 The form factors 

We will use the abbreviations for the five master integrals, used here and in the following for the 
s-channel contributions: 

Ao = Ao(m) (3.20) 

5o = 5o(0,0;s) (3.21) 

Bt = Boim,m;t) (3.22) 

Co = Co(m, 0, m; m^, m^, s) (3.23) 

Ci = Co(0,m,0;m^m^s) (3.24) 

Dq = Do("^) 0, m, 0; m^, m^, m^, m^, t, s) (3.25) 

together with the function 

C4 = Co(m, 0, m; m^, m^, if:). (3.26) 



The latter may be expressed by Aq and Bt, see (IA.8I) . This function contains the infra-red sin- 
gularities and we decided to keep it explicitely as it is also done in Loop Tools. Further we 
introduce 

w = (3.27) 

s 

^ = YZl^' (3.28) 

s 

y 1 4m.^ ■ (3.29) 

t 

In terms of (13.201) - (13.261) the results for the amplitudes are given in the following. They may 
also be obtained in FORM format from [43]. The explicit expressions for the master integrals are 
discussed in Section|4| The form factors from self-energies and vertices are: 

(3.30) 



F;'" = 2 - xz] + B.i|i(l + 2) + <i - 4] 

\ s a — 3 s 
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2 2 4rr7 
F^* = Ao^x[—- -{d- 4)] + 5,— a;[l - {d - 4)]. (3.32) 
ms^ a — 3 

For the box diagrams: 

= (Ff « + , (3.33) 

and the six independent box form factors for that are: 

— v\. A^v{v^ + 2w + vwx) + 2zw^'\}\ \Bt{w^y + v)z + i?o'i^w(l + ra)] 

C4 



d-3 



[(1 - Zz)yz^ - v{{\ - 6z)yz - 2(2 - 3^)^) -v\A+ {3y - 10)z) - Av^ 



+Civ[-^{2 -z-z'^ + v(2w+ + z)} - 2(1 -z)- 2v{l - wx)] 
d — 3 

sv I 1 

-Do— < [ -2vv_v+ +z{2 - z- z^) + v{l + 4:z)z - 5v^z] 

2 a — 3 

+2^;(1 -3z + 2z'^) + v'^{A - bz)+ 2v'^ - (1 - zfz] ) , (3.34) 



pm ^ —fAd[w + 2v{l + wx)+w^y]--[Bt{2v + w^y) + Bo{l + 2vx)w] 
s \ s 



d-3 

+Ci|^^[(l - 2^;)(1 -z)- 2v'^\ - (1 + 2vxz){\ - z) 

^-2v{v{l -xz)-z)\- i)o^{^[(l -z){\^w- 2vz) 

-2v'^{w + z)]+v{l + 2v{w+ + z)- 2-2^) }) , (3.35) 

F3™ = ^(-5^§^». + ;^{C.-Z,4.4). (3.36) 

s \ ^ s 

y - o (^- + ^) - (y- - - o)] - 



+ v{2v — z) — z 



d-3' 2'^ ' ^ 2" 2 d-3 

- (1 - z)] + Do^[v -^(w + w+-z)]), (3.37) 
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Fg^^ = ^i^Ad[v + {l + vx)w + {l + 2w)y]- (3.38) 



2vm 

s 
2 

-[Bt{v + (1 + 2w)y) + Bo{l + vx)w] - 

-w+{y- - 2wy - z) - Ciw\- — - - (1 + vx+)\ + -Do^[- 



rf- 3 ^ ' ^rf- 3 ' 2 ^c^- 3 

Fg^^ = 1 Arfb + 2w?/ + vwx\ \Bt{v + 2wy) + Bqvwx\ 

s \ s 

^ -{v — z)(2t> + to?/) — 



d-3 



1 SIJ 1 \ 

+ + ^+ - 2^)] - ^0 yi^l^ - ^) - 1]) > (3-39) 

where we have further introduced 

Aold-2 

Ad = 2-:r^^ (3.40) 

s a — 3 

x± = l±x, (3.41) 
?/± = l±y, (3.42) 
t;± = l±v, (3.43) 
w± = 1±U7. (3.44) 
The small mass limit is easily obtained by putting z = 0, x = y = 1. 

3.2 Counter term contributions 

In this section we focus on the contributions originating from renormalization: the charge counter 
term, the mass counter term and the wave function renormalization are given in arbitrary dimen- 
sion. Not only their 1 /e and constant terms are needed in order to render the amplitudes from the 
diagrams in Fig. (13.11) finite, but also 0(e) terms combine with divergent parts of the unrenor- 
malized amplitudes to give additional finite contributions. Similarly, of course, in two-loop order 
the 0(e) contributions of the diagrams combine with the 1/e terms of the counter terms to give 
finite contributions. 

First we consider the charge counter term. Each diagram of Figure (13.11) has at its vertices a 
factor e, the electric charge. Renormalization in two-loops requires e to be replaced by e{l+5e/e, 
with the charge counter term 

Se ^ d-2Ao{m) 

e (47r)^/2 3 ' ^ ^ ^ 

While the introduction of the charge counter term results only in an overall factor, the introduc- 
tion of the mass counter term is more complicated. Every internal electron propagator in Fig. 13. II 
has to be replaced by: 

1 1 / 2iv?5m/m , 

' ^ ' ' ' (3.46) 
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with 



6m 
m 



(rf- l)(d-2) Ao(m2 



(3.47) 



the electron mass counter-term. This means that additional amplitudes F^^'''^ are obtained from 
the one-loop diagrams Dk: All fermion propagators are replaced according to (13.461) . but the 
higher powers of Sm/m are dropped. The contributions from the first powers of 5m /m lead 
to 'dotted propagators' with squared numerators. The second recursion relation given in (IA.7I) 
reduces the resulting Feynman integrals with dotted lines to master integrals. 

Since the mass renormalization conterterm 5m /m contains a 1/e pole, the one-loop master 
integrals Aq and Bq etc., resulting from the diagrams, are needed to order (9(e) in order to ensure 
that all finite terms are taken into account properly. The relations between the nine amplitudes 
of a given dotted diagram fulfil similar relations as those for the undotted diagrams. The only 
differences are the following: 



m 



+ F. 



D8,r 



D9,r 



m 



-pD9,r 
7 



-Qdm 



^D9,r 



3 ' 



2m F: 



+ 2dmFr' + 2dF' 



iD9,r 



(3.48) 
(3.49) 
(3.50) 
(3.51) 



All other relations remain unchanged. 

Now we present the contributions to the amplitudes for the dotted diagrams: 



^1 



self,r 



Fi 



vert,r 



Ao — {d-2)xz 



'd~3 



Bs-z[l - {d-3)x] , 
+ 6^)-(^-4)(2x 



(3.52) 



1 + (rf-4)) + 7-4x 



m^s 



[1 + 3-^] - Bs-xz[{d - 4)(2x + (rf - 4)) + 2x\ 

(X s 



-ivert.r 



{d-3)s 
Ao 4 

2^ 

m 



xzis 



4m^)Cc 



(3.53) 



^ + (rf-4)(l + (rf-4)(l-|)) 



3 + 2x 



4m 
- B,—x 



[{d -A){x - (d- A)xz) + 2xz] 



4m 



{d-3)s' 



+ x)(s -4m^)C4 



(3.54) 



For the box diagrams: 



-^box,r 



[VW 



rpD8,r , rpD9,r 
j ' j 



(3.55) 



and the six independent dotted box form factors for that are: 

pD8,r ^ ^-I(d-A)\2v^(v+(l + v+)-xvw + 2wz)-zv(v(2 + v+)-wy(l + 2w) 
m"' s I L V 
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2 4 

z y w+ 



+ 



d-5 



vz 



V (^+(1 + — vw + 2u'2;) + —{v{2w — Zz) 



v'^z'^ 



2v'^[v{4: + 2w+ + 3;^)+ 6(w+ - w) - 2xvw] 



+z{v + y)^ 
+2y{wl-w)) + 

+y Y"+ - ^v{l + i;+ + — ) + y(3 + llw + 12w2)J + —[v{%v - 2 - 4yw_) 
-V(2 + «;)]}+ -^o ^,^^^_3^ {^[?/- + ^(1 + %) + 2^;^ - (1 - 2w+|/)^] 
+2wy^z^ - Bt-^ |4(d - 4)[t;(T; + yw+) + y + + + 4t;(l - 

+5yw_ - ^) + ^ (^v + 2v^ - 2y{vw -?/))}- Bq ^^^^ ^ ^2v{w + w+) 



d-3 



+2xvw{l + 2v- z) + 2y{w^ + wl)^ + C4 

+ + 9y) - (5 - 2?/ + 2«;(2 - y))z]^ + Y ) ~ + "^^^ + ~ ^^'^^ 



— + (1 - 2«;+y)— J -«^ry I + Ciy |(o? - 4) - ^;(1 + 2w + 8w' - Av') 

+4:xv^w + 2y{w'^ + w+) - z{l - 4vw+ + 2v^ + 2y{w^ + w+) ) - 2^^vj + 

— 2i>(u' + 2vv+ — 2xvw) + z(^w+ + {3v — z){v + v+f^ | + Dovzsl^{d — 4) 

2 

vvl - I (^v(2 - 2w + 3v) - y{w^ + w^)) + y - t; (^2 + (2v + + 
+1 (^;(3 + 2«;) - y{w' + <)) }, (3.56) 



^0 1 



I (d - 4) 2v (2v{2 + wa;_ + ^) + 4w+y - wy_^ - 2z (2v{v + y) 
+(3 + bw) y^' 



vwy- + 3 — ^[f + w+y\ + 2v ov[6 + v) 



d-3 



rf- 5' 



-(3 - lOv + 4vx - 15y)w + 12y - 2z 4v^ - 2vwy + (6 + 9w)y^j | 

+^ %3(/_3) i^^' + + + 2y)} - {^(c^ - 4) [2^;y + (3 + 5w)y^'\ 

+vwy^ — z (2v{v + yw+) — (3 + 4w)y^^ | 

+5o^|(rf - A)[y^ - z{l + 2xvw)] + y- - z{l + 2xvw)^ 
2 ( 

~^^ s{d - 3) Y'^y-'^^\^^^ ~ yy-) + ^^(^ + ^(^ + ^+))] 

+z^[wy'' + 2vy{2 + w) + 2v'^] | 
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2v 

-Ci — {{d - A)[y_ - z{l - 2v{w+ - xw) - y)] + 2zv{l - wx+)] 
s 

+DoVz{{d - 4)[2vw+ + y]-y+ + z}, 



(3.57) 



^ d-2 2 „1 ^ 1 , . 

— — — T^—v^y ~ Bt-vwyz — Ci-^^ — -;^vw[v — z)yz + C\-v wz 



'd-3s^ 

Dq-sv'^wz, 



2(d-3) 



(3.58) 



-^{y-- z) + {v - -y )yz 



-vwy+ + v{v + y) - -W-y^z 
{w + ^){l + 2y) + ^y + 3v-^ 



w? s^Xd — 
+ (d-4)[l. 

\{w-2)y'z] 

+Bt^ I (v{-y^ + (1 + 2y)z) - yh) - {d - 4)t/S(l + v+-z)j 



-Bo{d - 3)- {y_ + 2vx{v+ - z){l - z) - 2v{y + v - z) - z{l - 2y)} 



s2{d-3) 



{{v{l + 2z) - wy^z)y - (1 - z)v} {v - z) 



—Ci— \vz — {d — A) {w -\ — )7/(l — z) -\ — {v — w) — vz 
+D,vz -2z)+^-^ + \vy^ -{d- 4)[|(1 - 2z) - |> , (3.59) 



— — |(o? — 4) —2v(v + w + xvw + ywj^)-\- 2zy' 
2v\l — 3v — 2w — 2xvw — y{A + 3w) + 



w 



{yw+ + v)+ 



d-3 



d-5 
2z{v - 2y^w)} + 



4m 



Bt—{vy^ + z{{d-3)y^w-v)} 

s 

Amv z 
+Bo^{d -3)[v-w + y- 2vwx{l - -)] + 

^ 2vm 1 , , . X ^ 2-?7m f . , , , ^ . 
O4 ; — -{v — z)[y- — z) + Gi <{d — A) Zvw — v + w — y + z[v + y+ 



s d — 3 

xvw) + zxvw\ — Do2vm < -{d — 4:){v + y) + vv+w+ -[y- 



vw — z 



(3.60) 



— — I (0? — 4) 2v(v — xvw + y) — 2zy 
m y L 



W- 



+ 



d-5 



[v{v^yw^)\-^ 
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2v (^v — 2xvw + 3y^ + 2zy (2v — y{2 — w)^ | {vw + z{wy— 

2v)} - Bt— {{d - A)zy^w^ + yvw - zy{2v - y)} - Bq— {{d - 'i){v{v + y) 
w + zxv'^w)} — Ci—— — — {^—wyz^ + v{w{v — y-z) — 2{v — z)z)^ + 

2vTfl 

Ci {{d — 4)[(t> + y)w — z{y — w{yx — y))] + vw — z{2v — xvw)} + 

Dovm {{d — A)[z{v + yw)] + vw{v + w) + z{v{v- — 2w) — yw) + z'^v} . (3.61) 
Finally we investigate wave function renormalization for the electron self-energy: 

= A{p^) + _ m). (3.62) 

The wave function renormalization is given by: 

dA 

and the 'undotted' one then reads 



B + 2m— |p2=^2 = l + 5Z, (3.63) 



5Z = { + Am-DB,{^, m, m^) ) (3.64) 



with 

dp 



DBo{0,m,m'^) = ^^\p2^^2. (3.65) 



The first part in (13.641) contains the UV divergence and the second the infrared divergence. Ex- 
plicitely (13.651) reads 



{d — 3) 4m^ 

1 r(3-f) 

(m2)3-| (rf-3)(rf-4) 
Co(m, 0, m; m^, m^, 0) 
(rf-3) ■ 



(3.66) 



The UV divergent part of the wave function renormalization cancels the UV divergence of the 
vertex, and a remaining IR-singularity will be compensated by soft photon radiation. It is worth 
mentioning that due to (13.661) we can also write 

5Z = (3.67) 

m 
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We now discuss the dotted diagrams. They are UV finite and the divergent contributions 
to the 'dotted' 6Z come only from the IR divergence. Therefore we write the wave function 
renormalization from the dotted self-energy 5Z'^ in terms of DBq (or Co, respectively): 



6 dfji 

5Z' = -2 — {d - 2) [6 + 4{d - 4) - (rf - 4)^] DBo{0, m, m^). (3.68) 



The true one-loop form factors contribute to the interference with Bom (as shown here ex- 
plicitely) as well as to the squared one-loop correction (not shown explicitely), while the dotted 
form factors contribute only to the former. 

4 The master integrals 

The five master integrals of massive Bhabha scattering are shown in Figure |4~T1 We collect here 
expressions for them valid in d dimensions, but also the necessary e-expansions. 




The resulting form factors are: 



J Q 



(3.69) 




Tlllm 



SE212m 



SE210m 




V311m 



B412m 



Figure 4.1: The five one-loop Mis. External solid (dashed) lines describe on- (off-)shell mo- 
menta. 
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4.1 One-point function 

The simplest master integral is the tadpole: ^ 



Ao{m) 



dfik 



d/2 j J^2 



V{l-d/2){m^) — 

- + (1-LJ + |(2 + C2-2L„ + L^) 



—m 

J2 



(6 + 3(2 - 2(3 - 3(2 + (2)L„ + 3L^ - L^) 



+ ... 



with the abbreviation 

= In(m^) 



(4.1) 



(4.2) 



Often, shorthand notations with m = 1 are used, and our tadpole formula then agrees with 
Tlllm as it is given in the Mathematica file MastersBhabha .m located at [44]: 



-Ao(l) = Tlllm 



= 1+ 



e+ 1 + 



(2 (s 



+ . . . 



(4.3) 



4.2 Two-point functions 

The two-point functions are 



So(m,M;p2) = 



IT 



d/2 



m2)[(A;+p)2-M2]' 



(4.4) 



There are two of them, -Bo(0, 0; p^) (coming from the reduction of box diagrams) and BQ{m, m; p^). 
In d dimensions, they have been determined in [45] and in [46], correspondingly: 



So(0,0;p2) = 



e^^-v^ r(2-f)r(f-i) 



,2^2-1) 



2^-3 r (^i^) 



Bo{m,m;p') = e^^-(m^)-(2-t) r ( 2 - | ) 



d 



1,2 



3 

2 ' 



d . ^2 

2 ' ^ 



Am? 



(4.5) 
(4.6) 



The e-expansion for i?o(0, 0;p^) is trivial, 

1 



So(0,0;p^) 



+ 2 - In(-p') + e 



4-|-21n(V) + ^ln'(-p') 



+ ..., (4.7) 



^We omit here and in the following the conventional scale factor (47r/i^)^; the scale factor would make the 
arguments of logarithms dimensionless. 
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and the one for BQ{m,m;p'^) may be determined by using a relation for contiguous hypergeo- 
metric functions 



1,2- f; 



3 . 

2 ' 



4m^ 



l-2e 



1 - 2e 1 



Aw? 



1,1 + e; p 

3 . 

2 ' 



Am? 



(4.8) 



and then expanding the transformed hypergeometric function [46-49]. ^ The result is: 



Bo{m,m;p'^) 



- + 2- Lm + rln{x) + e 



A + ^-2Lm + lLl+r 



(2\n{x) 



+e 



- ln{x)Lm - 2 ln(x) ln(l + x) + - ln\x) - C2 - 2Li2(-x)) 

+ Ll + C(2) - ^C{2)Lm + r (-2C(2) + C(2)L„ 

-2C(3) - 4Li2(-,T) + 2Li2(-x)L„, - 2Li3(-.T) + AS,^2{-x) 
+2 ln(l + x){C{2) + 2Li2(-a;)) + 2 ln(a;) ln^(l + x) 

-(4 - 2L„) ln(x) ln(l + x) - ln(x)^ ln(l + ^) + ^ + 



;i - \n{xr + (4 - 2L„ + hi - ic(2)) ln(a;)) 



+ 



(4.9) 



with 



X 



r — 



1-^ + 1 



— 4m^ ' 
1 + x _ 1 
1 -X ~ b' 



1 + 6' 



(4.10) 

(4.11) 
(4.12) 



The e-expansions may also be determined by the method of differential equations [50, 51] and 
are then naturally expressed in terms of Harmonic Poly logarithms [44,52]. With m = 1 we 
have [44]: 

So(0,0;p^) = SE210m(x) 

= - + 2 + i/[0,x] + 2//[l,x] +6|(4-C2/2 + 2i/[0,x]+4//[l,x] 



+H[0, 0,x] + 2 {H[Q, 1, x] + H[l, 0, x]) + 4 H[l, 1, x]) } + . . . (4.13) 



So(m,m;p^) = SE212m(x) 



(4.14) 



- + 2 + H[0, x] + e|((-8 + 8 X + C2 

e 1 — X I 



^We thank M. Kalmykov for the FORM code hYpergeometric2Fl for an automatized derivation of the e- 
expansion. 
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+3 X C2 - 4 (1 + x) H[0, x] + 4 (1 + x) H[-l, 0, x] 
-2 H[0, 0, x] - 2 X H[0, 0, x])/(2 (-1 + x))) } 

+e2|(((l + x) (8 - 16/(1 + x) + 3 C2 - (2 C2)/(l + x) 
+ (5 C3)/3 + (2 C3)/(3 (1 + x)) - 2 C2 H[-l, x] 

+ ((-8 + C2) H[0, x])/2 + 4 H[-l, 0, x] - 2 H[0, 0, x] - 4 iy[-l, -1, 0, x] 
+2i7[-l,0,0,x] +2i7[0,-l,0,x] -/7[0,0,0,x]))/(-l + x))} + ... 

With the Mathematica file HPL4 .m., also located at [44], the corresponding expressions in 
terms of poly logarithms may be derived from (14.131) and (14.141) . 



4.3 Three-point functions 

There are two three-point functions, Co{0,m,0,m^,m?,p'^) and Co{m,0,m,m^,m^,p'^), with 
the definition (p = pi + P2)'- 



Co(mi, m2, ms, m\ m^,p^^ 



(4.15) 



[/c^ — mf] [{k + — ml] [{k + pi + P2Y — rnf^ 



As shown in (IA.8I) Co(m, 0, m, m^, rn^.p^) is not a master integral. The UV-divergences of Aq 
and i?o in (IA.8I) cancel, and the factor l/((i — 4) represents the IR-divergence of this vertex 
function."^ Due to the additional factor of 1/e, we need Aq and Bq up to O(e^) for a Co of order 
e. As discussed above, a separate control of IR divergences is often quite helpful in applications; 
therefore the explicit use of Co(m, 0, m; m^, m^, s) is recommended and we reproduce it here 
for completeness (see also Equation 39 of [46]): 



Co(m, 0,m, m ,m ,p ) 



ln(x) 



ln(x)[ln(— p^) + ln(— x)] 
1 



(4.16) 



2Li2(-x) 

2^ 



+2Li2(-x) - - In^(x) + C2j + e(^- ln"(a;) + 

ln2(x) + C2 [ln(-/) + In(-x)] + ^ ln(x) [ln(-p' 

+ In(-x)]' - I ln(x) + 45i,2(-x) - 2Li3(-x) - 2(3) + ■ ■ ■ } 

The vertex master integral Co(0, m, 0, m^, m^, p^) md dimensions is finite for small e; it has 
been derived in [54]: 



Co(0, 171,0,171 ,m ,p ) 



-e^^-(m2)f-3r(2-^ 



1 



2(rf-3) 



2F1 



1,1 



' ' 1 - 



p^ 



; 4m^ 



*The loop functions ^0 and Cq used here deviate by an overall sign from conventions of e.g. [39, 53]. 
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(4.17) 

Concerning the expansion with respect to e, the two coefficients of the 2-^1 -functions depend on 
Lm and ln(— p^), respectively. By eliminating Lm according to 

Lm = ln(-j92) - (21n(l -x) -ln(a;)) (4.18) 

we obtain 

Co(0,m,0,m^m^p2) = _J_|i 1^2(3;) + 2Li2(a;) + 4(2 - e ln(a;)=' - ^ ln(x)2 

+ [ln(V) - 2 ln(l + a;)][i \^{x) + 2Li2(x) + 4(2] 

+ ln(x)[2Li2(x) -2Li2(-x) + 5C2] -2Si,2(x2) +4Si,2(-x) 

+8Si,2(x) + 2Li3(-x) - 5C3) } + . . . (4.19) 

In terms of HPLs, the function reads for m = 1 [44]: 

Co(0,m, 0,m^m^p2) = - V311m[x] 

(4C2 + i^[0,0,x] + 2/7[0,l,x]) 



(1 -X2) 



e X 



5 C3 + 8 C2 x\ - C2 ^[0, x] + 8 C2 x\ 



(l-a;2) 

+2 E\-\, 0, 0, x] + 4 E\-\, 0, 1, x\ + i/[0, 0, 0, x\ 
+2 [0, 0, 1, x] + 2 if[0, 1, 0, x] + 4 if [0, 1, 1, x] 

+2/f[l,0,0,x] +4fr[l,0,l,x] 



(4.20) 



4.4 Four-point function 



In Loop Tools notation [39], the four-point master integral in d dimensions with two photons 
in the s-channel is : 

Box(t, s) = D^ym ,m ,m ,m ,t,s,m ,0,m ,0) 

= - / — (4 21) 

We first give the e-expansion obtained from a representation based on generalized hypergeomet- 
ric functions; see Subsection 14.4.11 Here we collect and complement results presented in [54] 
and [55]. Given the general result for the box diagram in d dimensions, the coefficients of 
the e-expansion are naturally obtained in terms of one-dimensional integrals. Alternatively we 
consider in S ub section 14 . 4 . 21 the method of differential equations, which also yields the coeffi- 
cients in terms of one-dimensional integrals. These can, however, systematically be presented 



20 



in the form of generalized harmonic poly logarithms, which makes this form quite attractive if 
one prefers 'analytic' results. Finally, in Subsection 14.4.31 we add a representation in terms of 
a two-fold Mellin-Barnes integral, which appears to be quite elegant and has the advantage that 
the integrand is free of singularities even in the physical domain. 



4.4.1 Hypergeometric functions 

A closed expression for the box function valid in d dimensions is known from [54]. In this case 
a first order difference equation with respect to the dimension d was solved. ^ Other difference 
equations use as parameter the powers of the propagators , see e.g. [56-58]. The general result 
of [54] reads: 

T. / N 4m'^"^ ^/ d\^fd-3_,3d-2 t 
e-'^^Box(t, s) = -r 2 - - Fa , 1, 1, -, ; z 



4m'^ * „ / d 



(c/-3)s(t-4m2) V 2 

0F(_,)^ r(^)r(2-f) ^ / d-3 ^ l_ d-l 

2<^-^m4~s (t-4m2)r(^^) 2 ' '2' 2 ' t-4m2'^ 4m2 



4m2 



d-l d-2 ^1 

~2~- ~' S 



u ^ s 



(4.22) 



with 



4m2 u ^ 
s t — Am? 

The two photons are in the s-channel. Naturally the cuts of the diagram are different for the 
t-channel case, which means that the hypergeometric functions are to be evaluated in different 
domains of analyticity. In (14. 221) . e.g., the imaginary part of the diagram comes only from the 
coefficient {-sY'^'^^''^ of F^. 

The Appell hypergeometric functions in terms of their integral representations are: 

I d-l \ d-3 dtt'^ 



'd-3 



3 d-2 



2 '^'^'2' 2 '"^'^ 



and the Kampe de Feriet function is: 



Ut 



^/T^{l-tx)- 



— 2 



d-3 d-2 



1 - tx' 



pl;2;l 
i;i;0 



2 ■ 2 
d-l . d-2 
L 2 ■ 2 



1; 1; 



d-3 dtt^ , d-3 d-2 

^Fi{l,^^,^^,zt). 



(4.25) 



(4.26) 



lol-ty^'^' 2 2 

^We just mention that in [55] also a Feynman parameter representation for Dq was given, including terms pro- 
portional to e. 
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See [59] for (14.241) . [55] for (14.261) . and (14.251) is obtained from the double integral representation 
of the F2-function [60]). With these representations we can derive the neeeded e expansion. Due 
to r(2 — |) in the prefactors of (14.221) . their e-expansion has to be done up to order e^. This 
can be performed by expanding the integrands. The numerical evaluation of the one-dimensional 
integrals of the e-terms works quite nicely in general. Nevertheless partial analytic results can 
also be obtained, see e.g. (14.271) and (14.661) . Based on [55] we also give an expansion of 
the integrals for the limit of small masses, i.e. —t ^ 4m^ (neglecting terms of O(m^) and 
0{mHn{m^) ). 

For the Fi -function the e-expansion is easy except for the analytic integration following the 
expansion. In [54] the analytic integration has been performed for an Fi -function in which one 
of the arguments is 0(e) and in [55] the corresponding transformation to obtain such a form has 
been described in detail. For the real part of Fi we thus have: 



'd-3 1 d-1 



x,y 



Re\ ln(S) 



m d — 3 

-e(Li2(l - AB) + Lis " f ) " 2Li2(l - B) 



+2Li3 
+2Li3 
+2 



A-B 

Ml - B) 

1-AB 
1- B 



A-B 



-2Li3 
2Li3 



1-AB 

Ml -B Y 

A-B 
1-B 



+Li2 

+ 
1 



Li ( MA-B 
[T^AB 
1-B 



Li 



1-AB^ 

A{1 - B] 



A-B 



^ln^(A)-C(2) 



1-AB 
Li2(-A)J ln(A) 
B-A 



In 



6 \1-AB 



1 W.N, 2 f B-A 

2 ^ ^ \1-AB 



1-AB 



(4.27) 



with 



A 
B 



x{s) = 
-x{t) 



m 



a + 1 
b-1 



m 



< 0, 



<0, 



(4.28) 
(4.29) 
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S — 4:1V? 

s 

t 



(4.30) 
(4.31) 



t - Am? 
Abbreviating (14.271) as (6 ~ 1) 

^1 1, \ ^, y) = -^(^ - 3) K + ^Fl + e'Fl] , (4.32) 

we obtain from (14.271) in the limit of small masses with r = — t/s, < r < 1: 





-ln(4)' 




-i). 

r 


Fl = 




-2C(2)-Li2(l 


Fl = 




-2C(2)ln(4) 


_Li3(l-i)-2C(3) 
r 



(4.33) 

For the F2- and Kampe de Feriet functions the same hypergeometric function 2F1 needs to be 
expanded: ^ 

2Fi(l, J-e,l-e,^i;) = J {l - 2e ln(l + v) + 2^ + v) + U2{-v)\ ] 

^ {(1 + v)~^' + 2^U2{-v)] + 0(6^), (4.34) 



with 



. = ^"V^ . (4.35) 

For the F2 -function we have to use 

z Av 

w = = (4.36) 

1 — xt (1 + vj^ 

and correspondingly for the Kampe de Feriet function 

Av 

w = zt= (4.37) 

(1 + vy 



and further in the integral (14.251) 

(4.38) 



+ ,z 



\-tx)-^ Av' 



^ Again a FORM code for the automatized derivation of the e-expansion by M. Kalmykov has been used. 
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and in the integral (I4.26I) 



-2e 



(4.39) 



It appears natural to introduce v as integration variable. But a more precise numerical integration 
results from an elimination of the singularity at t = 1 in (14.251) by the transformation 1 —t 



u . 



We then have in the considered order for (14.251) : 

'd-3 3 d-2 \ du 



2 '^'^'2' 2 



2 , 4m2 



with 



(4.40) 



z 

Av 



X 

4 



— 



t 



Am? / „ Am? , Am? . 

+ a/m2 + 1 



t 



< 1, 



(4.41) 



and 



U2 - ^ 



4m£' 

t ■ 



U 



2 'im? 



> 0. 



— + + —[I —) 



Am? " 



For the following we write 

'd-3 3 d-2 

' '2'^"'"^'^ 



(4.42) 



(4.43) 



where is obained as 



:ln 



(4.44) 



and the higher orders must be calculated from (14.401) numerically. In the limit of small electron 
mass they are: 



^2 =-iln'(i) + C(2)-Li2(l--), 



m^ 



r 



6 m^ 



m^ 



Fi = \ ln^(:^) - ln(-^) f C(2) - Li2(l - \)\ + ^C(3) + S,,,{1 - -J - U,{1 - 



(4.45) 



Similarly we perform the calculation for the Kampe de Feriet function (14.261) : 



-^1:1 



1;2;1 




2 • 2 

d-1 . d-2 
2 • 2 



, 1; 1; 



z,y 



d-3 



2 
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du 



with V = vqu'^ and 



1 1 

+ 



1 + biu 1 — biu 



hi 



1 1 

+ 



1 + h2U 1 — 62W 



(4.46) 



^^0 



hi 



1 + 



z z 
4 



^1 



1 + Jl- 



— <^ 1 and 

Vl 



^^2 = — , 
Vl 



— = VvqVi- 

V2 



As above for the F2, we formally write for the Kampe de Feriet function: 



pl;2;l 
i;i;0 

where JC^ is obtained as 

;i-&i)(i + &2 



1; 1; 



d-3 . d-3 

d-1 . d-2 . _. ^; ?/ 

2 ■ 2 ' ' 



ci-3 



+ eK^ + e^K^ + ■ ■ ■ ) 



In 



'l+hi){l-h2) 

Again, investigating the small mass approximation, we have 

= 3C(2), 
= 7C(3). 



(4.47) 



(4.48) 



(4.49) 



(4.50) 



Finally we see that the expansion in e of the F2- and Kampe de Feriet functions becomes easy 
with the representations (14.401) and (14.461) . To sum up our results, we have 



Box(t, s) 



2(m 



2\-e 





— 4m2) 






s{t 


- Am?) 






s{t 


— Am?) 



r(e) [F!^ + eFl + e'Fi + ...] 

r(l) 



re 



r (1 - e) [Re [F^ + tFl + e^F^ 



■}]• 
(4.51) 



As we see, in the limit of small electron mass the 1/e-terms of the F2- and Kampe de Feriet 
functions cancel. 

It is a very appealing fact to have the closed form of the box function as an analytical expres- 
sion in d dimensions. So far, however, only partial analytic results were obtained for the terms of 
order e , but as we observe already from (14.271) . the results become quite lengthy if one prefers to 
present them in this form. Beyond that simple expressions in the small mass limit were obtained. 
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4.4.2 Harmonic polylogarithms 



An alternative approach in terms of solving a differential equation for the box [61] yields in a 
natural manner harmonic polylogarithms. For the purpose of checking (in particular also numer- 
ically) and comparing, we repeated the calculation of [61] and shortly sketch the procedure. 

To be explicit, we consider the Bhabha box diagram with two photons in the t-channel, as 
in [61], and the electron mass being set to 1; the analytical continuation to the s-channel is 
evident here. One may derive the differential operator 



d 



V1+V2 



d 



ds 2 ' ' s + t-A'""' dv% 
which applied to the one-loop box yields a differential equation: 

dBox(s,t) 1 



(4.52) 



ds 



2(-4 + s)2st(-4 + s + t) 

-4 + s)t{-2s^ + 4(-4 + t) + s(12 + (-6 + d)t)) Box(x, y) 

-2(-4 + s) (-4 + d) (-4 + t)t V311m(?/) 

+4st(-3 + d) SE212m(x) 

-4(-3 + rf)(-4 + s)(-4 + s + t) SE210m(?/) 



-2(-2 + d)st Tlllm 



(4.53) 



where 



X 



^1^477-1 

v/r^47t-i 



or 



(1 



X] 



X 



yy 



y 

The subdiagrams Tlllm, SE212m, SE210m, V311m are given in the preceding sections. 
Expanding now the differential equation (14.531) in e and introducing the ansatz 



(4.54) 
(4.55) 

(4.56) 
(4.57) 



Box = const B412m(a;, y) 

Bo + eBi + 



e 



(4.58) 
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we may iteratively solve a system of differential equations which differ only in the inhomoge- 
neous terms: 

, ^ = -r. :^B,{x,y) + C,{x,y). (4.59) 

ax x[l — x^) 

More details are described in the literature, e.g. in [61]. 
The result is 

2xyH(0,x) 2 , 
^-1 = n 2V1 ^2 = :r^ HiO,x), (4.60) 

where H{0, x) = ln(x) has been introduced, 

Bo = J—— HiO, x) (if (0, y) + 2H{1, y)) , (4.61) 

and finally 

Bi = Z^-— Ici--, 0, 0, x) + Gi-y, 0, 0, x) 

st^/i-A/s [ y 



-2 {g{-K -1, 0, x) + G{-y, -1, 0, x)) 

- {g{-^, 0, a:) + G{-y, 0, x) - 2H{-1, 0, x)) [H{0, y) + 2/7(1, y)] 

- x) - G{-y, x) + ff(0, x)) (0, 0, y) + 2/f (0, 1, y)] 

- {^G{~ x) - 3G(-t/, x) - ^i/(0, x) - 2i7(-l, x) - 4i7(0, y))C2 
-2 {h{1, y)H{0, 0, y) - HiO, y)H{0, 1, y)) 

-2 0, 0, x) - 2H{-1, -1, 0, x)) - 2H{0, x) [i/(l, 0, y) + 2H{1, 1, y)] 

+i/(0, 0, 0, y) + 2H{1, 0, 0, y) - 2(3 } • (4.62) 

The functions G are generalized harmonic poly logarithms [61,62]. For the calculation of Bi we 
used the relations 

G(-y, 0,0,1) = -C2//(0,y)+//(0,0,-l,t/)-i/(0,0,0,y), (4.63) 
G{-y, -1, 0, 1) + G{--, -1, 0, 1) = -^C3 - C2H{-l,y) + H{1, 0, 0, y) + H{0, 0, -1, y). 

y 2 

(4.64) 

There is a difference in the coefficient of the term H(0, x)(^2 w.r.t. [61] due to different choices 
of normalization, see also (14.601) . 
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In order to check the results, we evaluated both representations of Bi numerically (for the 
photons in the s-channel). For s = 10^, cos 6* = 0.4, m = 1, agreement to nine decimals was 
achieved: 

5i(s = 10^ cos^ = 0.4) = 4.43779985 10"^ - 1.61529999 10"^ i. (4.65) 

The most difficult part of the numerical evaluation of (14.621) is the calculation of G{—y, — 1, 0, x), 
in which case a principal value integral has to be performed with the above parameters. The 
imaginary part obtained from (14.621) therefore agrees only to 7 decimals with (14.651) . Following 
[55], a simple formula for the imaginary part of Bi can be derived, which is indeed simpler than 
what is obtained from (14.621) : 

$5(50 



[2U2{1 + xy) + 2Li2(l + y/x) + AU^^-y) + In'(-x) + — 
+21n(y)[ln(s) + 21n(l + ?/)]). (4.66) 



This yields the imaginary part of the above number. The numerical calculations were performed 
with Mathemat ica and Maple, respectively. 

4.4.3 Mellin-Barnes representation 

Finally, we derive a Mellin-Bames representation for the QED box integral, again with two 
photons in the s-channel. The Mellin-Bames representation reads for finite e: 

Box{t,s) = ^.^—^^-^^—^ dz, dz2 (4.67) 



r[-2e](-t)(2+^) (27ri)^ 



—too 



r'[-l - e - - Z2 



r[2 + e + ^1 + Z2]T^[1 + ;2i]r[-^i]r[-^2] 



2r r[-2-2e-2;zi] 



T[-2-2e-2;zi - 2z2] 

A derivation may be found e.g. in [63]. Starting from this Mellin-Barnes integral, one has to 
perform an analytic continuation in e from a domain where the integral is regular into the vicinity 
of the origin. The singularity structure near e ~ is obtained by means of the Mathemat ica 
package MB [64]. We obtain the result in terms of the following one- and two- dimensional 
integrals: 

II = -^TT- / dz,i — ] ^ ^ '\ (4.68) 



and 



-|+ioo 



12 = ^^ y_._.oo "^^^ W ' r[-^i]r[-2(i + zi)]r^[i + ^i] (4.69) 
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In terms of the conformally mapped variable 



II = ziT.T^My)- (4-71) 



y = \ , (4.70) 

^l-4mVt + l 

the first integral 1 1 in (14.681) can be performed analytically to yield the well known result 

i 2y 

im?s 1 — y 

The final result for the Box then reads: 

Box(t, s) = -ill + ln(-s)ll + e Q [C(2) - In^(-s)] II - 212^ . (4.72) 

The first two terms are in evident agreement with (14.601) and (14.611) . The double integral 12 
in (14.701) is not easily evaluated analytically, although we know the answer from (14.621) . The 
MB package yields fairly precise values in the Euclidean region (s < 0). In the Minkowskian 
domain (with s > and (— s)^^ = s^^exp{—i7izi)) our experience with Mathematica is that 
the built-in function NIntegrate with MaxRecursion — > 12 gives easily a precision of 
nine decimals. An alternative is the expansion at small m and fixed value of t. With 

2 

— m 

mt = (4.73) 

r = ^, (4.74) 

we have obtained a compact answer for 12 with the additional aid of XSUMMER [65]. The box 
contribution in this limit becomes: 

Bi = ^|4C3-9C2ln(m,) + ^ln='K) + 6C2ln(r)-ln2(mi)ln(r) (4.75) 

+ -ln^(r) - 6C2ln(l + r) + 21n(-r)ln(r)ln(l + r) - ln^(r)ln(l + r) 
3 



+21n(r)Li2(l + r) + 2Li3(-r)} + 0{mt). 



5 Summary 

A calculation of Bhabha scattering for the luminosity measurement at ILC is promoted by sev- 
eral groups, aiming at a precision of 0.01%. With this study, we provide a publicly available 
program for the one-loop electroweak Standard Model corrections. Further we collect all needed 
expressions for the factorizing one-loop QED corrections. They are necessary ingredients for the 
full two-loop calculation of Bhabha scattering. 
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A Reduction of tensor and scalar loop functions to master 
integrals 

We strictly apply here dimensional regularization, i.e. infrared as well as ultraviolet singularities 
are given in terms of only one pole in e = {A — d)/2. 

After using DIANA [36] and FORM 3.1 in order to express the Feynman diagrams in terms 
of tensor integrals, we have to express the latter ones by scalar integrals: writing the m}^ scalar 
amplitude of a diagram formally as 

n,l i,n,l 

where momenta pi^^ are the 'chords', i.e. momenta in the propagators ci = [k — pif — mf, 
with k the loop momentum. The generically denoted n-point scalar, vector, and tensor integrals 
In I, In I, In I wiU bc transformed into scalar integrals with shifted space-time dimension d, which 
are then reduced to scalar integrals in generic dimension by means of recursion relations [66-68] . 
The indices / in dA.ll) stand for 'dots' on lines /. The reduction to scalar integrals reads: 

n-1 

fJ- t[' 
n. n.; J I - - 



/a IL—-L 
KU^''^'''' = -J2p'^^^^nl 
r=l i=l 

/d " -r 

K n ^r^'^'-' = Yp' p) -^.^ - C ' (A.2) 



r=l *ii=l 

where [d+\ is an operator shifting the space-time dimension by two units , riij = (1+Sij)\, riiji = 
(1 + 6ij + Sii + 6ji - 6ij5u6ji)\ and 

/d " I f d'^k 

r=l ^ 

is the original scalar integral with additional powers (dots) of the i-th and j-th propagators. 
The case with no dots is formally obtained by putting j = I = 0. Having reduced the tensor 
integrals to scalar integrals, the generic space-time dimension d needs to be re-established and 
the dots to be removed. For this we use the recurrence relations first proposed in [67], which 
are complementary to those obtained via integration by parts [69,70], and later simplified and 
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extended to zero Gram determinants in [68]. With Yij 
determinant 



-{pi — pj)"^ + mf + m'j and the Cay ley 








1 


1 . 


. 1 
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Yn 


Yu ■ 


■ Yin 
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Yu 


Y22 ■ 


■ Y2n 
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Yin 


Y2n ■ 


Y 



(A.4) 



the so-called signed minors (^^^^ ') determinants where the rows ji,j2, ■■■ and columns 
ii, Z2, ••• are erased from the Cayley determinant Making successive use of the following 
three recurrence relations leads to scalar master integrals Aq, Bq, Cq and Dq in d dimensions: 



1=1 




(A.5) 



(A.6) 
(A.7) 



These relations are applied in a FORM program one after the other: the (IA.5I) reduces the dimen- 
sion and the index of the j^^ line, the (IA.6I) reduces the index of the j*^ line without changing the 
space-time dimension. The third relation (IA.7I) also reduces the space-time. The operators i+, j+ 
raise the power of the corresponding propagator by one unit, while reduces the power of the 
A;-th propagator by one unit. 

For Bhabha scattering in particular there is one subtlety: there occur zero Gram determi- 
nants and for this case special care must be taken. The occurrence of zero Gram determinants 
(e.g. On = 0) is discussed in [68]. Effectively a zero Gram determinant reflects the kinemat- 
ical boundaries of phase space where a given n-point function can be expressed through scalar 
integrals of lower rank. A typical example of such simplifications is 



Co(m, 0, m; m^, m^, s) 



4m2 



d-2 Ao(m2) 2d -3 
d- A m? d- A 



BQ{m, m; s) 



(A.8) 



(IA.8I) means that this Co is strictly speaking not a master integral. For practical reasons, however, 
we include it in the list of master integrals: see the discussion in Section |431 It is instructive to 
derive (IA.8I) in order to demonstrate how the procedure works. Setting the momenta of the 



■'Note here the additional overall sign (-i)ii+J2+-+n+i2+.. 
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incoming massive lines to pi and p2 (the third momentum q = — pi — P2, q'^ = s) and the 
integration momentum on the massless line (no. 3), then the chords are, respectively, — Pi,P2 
and 0. Correspondingly we have for the Cayley determinant 





1 
1 
1 



1 

-s + 2m^ 




1 1 

-s + 2m^ 

2m2 





Sis 



4m^ 



(A.9) 



Apparently ([J)^ = 0. Applying (IA.7I) with d ^ d — 2, we obtain 

(^-4) ()3/|=-E 



k=l 



{d-2) 



(A. 10) 



where 03 =©3 = and = -s{s - Am'). 

Now we have expressed our three point function already in terms of a two point function with 
two massive lines, however md — 2 dimensions, i.e. 



{d - = 



(A. 11) 



and we must increase the dimension again with the intention to obtain an integral with nonva- 
nishing Gram determinant. The relevant relation to be used is (29) in [68] 



-jid+2) ^ _ja 



(A.12) 



which in our case yields 



r(rf-2) 



-2 1+J, 



2 ! 



(A.13) 



i.e. a two point function in generic dimension with a dot on one of the two massive lines and 
we have to remove the dot from the line. In this case we have ()2 = —2s and (g)^ = —s{s — 
Am'), i.e. both Gram determinants are nonvanishing and we can apply (IA.6I) . which yields 
straightforwardly (IA.8I) . 
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